A METHOD AND APPARATUS FOR MAKING 
A CHANNEL ESTIMATE 


BACKGROUND OF THE INVENTION 

5 

1. Field of the Invention 

The present invention relates to channel estinnation in a telecomnnunication 
system, and more particularly, data-aided channel estimation. 

10 2. Description of Related Art 

It is very well known that coherent methods of channel estimation give 
about 3 dB of gain over non-coherent methods, provided that perfect knowledge 
of the channel fading parameters (amplitude and phase) is known. In practical 
applications, such perfect knowledge is never achieved, and a receiver must 

15 estimate the fading parameters. This estimation is prone to errors, and thus, the 
"promised" 3 dB gain may rapidly start to vanish as the accuracy of the channel- 
estimation algorithm decreases. This problem appears in communication systems 
where channel fading and coherent detection are present (e.g., CDMA2000, 
3GPP, and ARIB - where the channel parameters are usually estimated from a 

20 pilot signal). Common losses from channel estimation errors in these systems 
can range from 0.1 dB up to 3 dB, depending on the specifics of the channel and 
the system. 

Under consideration is the problem of estimating the complex-valued 
channel gain (also called "complex fading coefficient") experienced by a BPSK- 
25 modulated symbol transmitted over a single path of a (possibly multipath) fading 
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channel. A complex discrete-time received sequence is generated by 
demodulation (e.g., correlator or matched filter) and sampling of each path. For 
symbol / , the effect of fading on the received signal, , can be modeled as 

5 where is the complex-valued fading coefficient, e {±1 } is the transmitted 

data bit (which can be either +1 or -1), and is the complex background additive 

white Gaussian noise with zero mean and per-component variance . It is 
assumed that the fading is sufficiently slow so that the cannel gain is 
approximately constant over consecutive symbols. Hence, the symbol subscript 

10 on the channel gain a is removed. 

Conventional DS-CDMA channel estimation is performed through the use 
of a reference signal known as a pilot. The pilot may be transmitted in a number 
of different ways. One approach is to provide a channel, separate from the data 
channel, exclusively for the pilot signal. This method is used by both the 

1 5 European (UMTS) and the North American (CDMA 2000) third generation 

wireless systems. A second approach is to time-multiplex pilot symbols with data 
symbols. This approach is used, for example, in the Japanese third generation 
wireless system (ARIB). Although the above two methods have fundamental 
differences, the underlying concept is the same. 

20 With pilot-aided (PA) channel estimation, the pilot data bit, > *s known to 

the receiver. Without loss of generality, we assume that 1 for all pilot 
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symbols. Therefore, for pilot symbol / , the received signal statistic y^^ is 

5 where is the complex channel gain for the pilot signal, represents the 

background noise and is Gaussian, and p represents that a variable is based on 
pilot symbols. Typically, the transmit energy of the pilot signal is kept as small as 
possible in an effort to minimize the consumption of battery power and added 
interference. Hence, the pilot signal does not necessarily have the same energy 
10 as the data signal. Since the goal is to estimate a and not , the following 

weighting is applied to the channel gain of the pilot signal 

15 where y3 is a known, chosen design parameter. 

For each received pilot symbol, a simple individual estimated realization of 

the channel gain, a^ ^, can be formulated as 

^^.p ^yi,p • 

20 

Data-aided (DA) channel estimation offers an alternative approach by 

making use of the data symbols in addition to (or in lieu of) the pilot symbols. The 

difficulty with DA estimation is that the data symbols are not known a priori fas is 

the case with pilot symbols), which makes the estimation more challenging - and 
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more noisy too. DA channel estimation can be implemented in conjunction with 
PA channel estimation by means of a stage-by-stage iterative procedure. In the 
first stage, a PA-only channel estimate is generated. This estimate is then used to 
detect the data symbols, which in turn are used as new "pseudo-pilot" information 
5 to revise the PA channel estimate. This process can be repeated, and may be 
implemented as a loop within the decoding stage of the receiver. Typically, the 
information in each data symbol is estimated, then the information is removed and 
a channel gain estimate is generated using an averaging window. 

DA channel estimation is an efficient way to assist channel estimation 
10 because it makes use of data information that is already available at the receiver. 
The goal is to reduce the pilot symbol overhead and/or reduce the required 
transmit energy of pilot symbols. However, conventional DA channel estimate 
methods are computationally intensive and exhibit larger than desired variance. 


1 5 SUMMARY OF THE INVENTION 

In the method of the present invention, estimated individual realizations of 
the complex-valued fading coefficient, commonly called channel estimates, are 
generated. According to the inventive method, these realizations are easy to 
generate and can produce a channel estimate that exhibits a smaller variance in 
20 comparison to conventional methods. 

A level of confidence for each possible value of a transmitted data symbol 
is determined based on a received data symbol corresponding to the transmitted 
data symbol. The confidence level associated with a particular possible value is 
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the level of confidence that the transmitted symbol was the particular possible 
value. Using the confidence levels, a channel estimate is generated. 

Unlike conventional data-aided channel estimate methods, the method 
according to the present invention does not require making an explicit calculation 
(called a hard decision) of an estimate of the transmitted symbol. Instead the 
present invention offers a soft decision alternative that does not require 
performing any maximizations. Consequently, the methodology of the present 
invention offers an easy means of determining a data-aided channel estimate that 
exhibits a smaller variance in comparison to conventional methods. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention will become more fully understood from the detailed 
description given herein below and the accompanying drawings which are given 
by way of illustration only, wherein like reference numerals designate 
corresponding parts in the various drawings, and wherein: 

Fig. 1 illustrates a plot of the confidence function hiy,j) versus the log- 
likelihood ratio A, ;and 

Fig. 2 illustrates an apparatus for making a channel estimate according to 
one embodiment of the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In the method of the present invention, estimated individual realizations of 
the complex-valued fading coefficient, commonly called channel estimates, are 
generated. According to the inventive method, these realizations are easy to 


generate and can produce a channel estimate that exhibits a smaller variance in 
comparison to conventional methods. 

Using the received signal y,^^ , the channel estimate is defined as: 

^,4^Ky,Jy,,u- (1) 

where a, ^ is the channel estimate, h() is any predefined function that can 
be designed based on a specific constellation and channel being used, i 
represents the Ith estimate, and d represents that the parameter or variable 
pertains to data symbols (as opposed to a pilot symbols). 

For the purposes of discussion only, the method of the present invention 
will described for the bi-phase shift keying (BPSK) constellation and an over-the- 
alr communication channel. In BPSK, the transmitted symbol is either +1 or-1. 
However, from the following disclosure, it will be understood that application of the 
method of the present invention Is not limited to a particular constellation or 
channel. In a preferred embodiment for BPSK modulation, h{y, ^)\s defined as: 

h(y,,) = P(x, =+'i\y,,,)-Pix,=-l\y,,), (2) 

where P(x, = +l\y,^J is the a posteriori probability that a transmitted data symbol 
X, =+1 was transmitted conditioned on the observation of the received data 
symbol y,^^ corresponding to the transmitted data symbol x, , and where 

= -^\y,,ci) is the a posteriori probability that a transmitted data symbol x, = -1 
was transmitted conditioned on the observation of the received data symbol y,^ 


corresponding to the transmitted data symbol x, . Here, P{x^=+l\y, j) represents 
a confidence level, based on the received data symbol y, ^ , that the 
corresponding transmitted data symbol x, had a value of +1 . And, P{x, = -\\y^ ^) 
represents a confidence level, based on the received data symbol y,^, that the 
corresponding transmitted data symbol had a value of -1 . 

Using Bayes' rule, equation (2) can be rewritten as equation (3a) below: 


where p{y, j) represents a probability density function of the received statistic 
evaluated at >>, ^ . 

A similar expression exists for P(x, = -ll^,^) . 

ffe.-.|,„) = ^fa^l-^7'>f.°-'> (3b, 

Under the assumption that = +1) = />(x, = -1) = 0.5 using the Law of 
Total Probability, /?(j,,Jis given by equation (4) below: 


P{y,,a) = -iPiy,Ax, = +1) + p{y,^,\x, = -1)) (4) 
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The well-known log-likelihood ratio (LLR) is defined as: 


(5) 


where ln() represents the natural logarithm. It is known that y conditioned on xis 
a complex Gaussian random variable with mean ax and per-components variance 
(7^ , (i.e., noise), where a-- is determined according to any well-known technique. 
Therefore, the LLR is given by equation (6) below: 


^(y)^My)-^-,(y) = 


(y-a) 
2ct' 


2 A 


(y + a ) 


2 A 


2a* y 


(6) 


Combining Equations (1) - (6) results in: 


a, J = 


y 


i,d ' 


(7) 


Returning to Equation (2), h { y ^ ) js given by: 


(8) 


A plot of this function with respect to X is given in Figure 1 . The function is odd 
and is bounded by ±1 , and bears a strong resemblance to the signQ function. 

As will be appreciated, h( j,^ ) represents the confidence that the 
transmitted symbol is a particular value in view of the corresponding received 
symbol y,^^. Stated another way, ^{y,j) indicates the strength or degree of 
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confidence that the transmitted symbol is a particular value in view of the 
corresponding received symbol y^^. 

Next, the overall data-based estimate is found by averaging the individual 
realizations of the channel estimate over a weighted time window: 

where K^ j is a weighting constant, IN^ + 1 is the window over which the 
estimate is averaged, and A'^^ is a number of samples. 

Unlike conventional data-aided channel estimate methods, the method 
according to the present invention does not require making an explicit calculation 
(called a hard decision) of an estimate of the transmitted symbol. Instead the 
present invention offers a soft decision alternative that does not require 
performing any maximizations. A simple evaluation of the LLR, a well-known 
receiver calculation already made in most receivers, is used. Consequently, the 
methodology of the present invention offers an easy means of determining a data- 
aided channel estimate that exhibits a smaller variance in comparison to 
conventional methods. 

Once the appropriate pilot-aided (PA) and data-aided (DA) channel 
estimates and their variances are obtained, they can be combined in an optimal 
manner. In the present invention, optimality is defined as minimum variance in the 
final estimate. The PA channel estimate a^"^ can be determined according to any 


well-known technique, and therefore will not be described. The variances al and 
crj of the PA and DA channel estimates a^^^ and a^^^^can be determined 
according to any well-known statistical technique; and therefore will not be 
described. The final channel estimate, a is a linear combination of the PA and 
DA channel estimates. 


where w^and are non-negative constants. Assuming that E[a^^^] = Ela^"^^] = a , 
where E[] represents the average value, the added constraint that 


ensures that £[a]= a . 

Under the assumption that the PA and DA channel estimates are 
independent, the variance of the overall estimate is 


To minimize this variance subject to the constraint in equation (11) and 
>v^,w^ being non-negative, :=1-M^^is substituted in Equation (12). Then, 

equation (12) is differentiated with respect to w^, set equal to zero, and solved for 

w^. The result is 


(10) 


(11) 


Var(a) = w^al + w]crl . 


(12) 



(13) 


and 



(14) 


A check of the second derivative confirms that the solution is indeed a minimum. 

Accordingly, by substituting equations (13) and (14) into equation (10), the 
channel estimate can be calculated using the PA channel estimate, the variance 
of the PA channel estimate, the DA channel estimate and the variance of the DA 
5 channel estimate. 

An apparatus for implementing the above-described embodiment of the 
present invention will now be described with reference to Fig.2. As will be 
appreciated from the forgoing, the apparatus of Fig. 2 forms a part of a receiver. 
Because the other components of the receiver are well-known, applicants have 
1 0 not illustrated and will not describe these other components for the sake of brevity. 

As shown in Fig.2, a shift register 10 inputs the received symbols from a 
demodulator (not shown). As alluded to above, also not shown are the well- 
known components for determining the per-component noise or variance <y^ , the 
pilot-aided channel estimate a^, the variance al of the pilot-aided channel 
1 5 estimate and the variance of the data aided channel estimate. An LLR 

calculator 12 receives the received symbol y- from the shift register 10 and the 
square of the standard deviation, and calculates the LLR of the received symbol 

according to equation (6). As will be appreciated, equation (6) requires, during 
a first iteration, an initial channel estimate as an input variable. In a preferred 
20 embodiment, the channel estimate based on the pilot symbols is used as the initial 
channel estimate, and then each subsequent iteration uses the channel estimate 
determined based on the combined data aided and pilot-aided channel estimates 
as shown in Fig. 2. 


Next, a confidence factor generator 14 generates a confidence factor 
h(y) according to equation (8) using the output of the LLR calculator 12. A 

nnultiplier 16 multiplies the confidence factor with the received symbol on a per 

component basis to obtain an individual realization of the channel estimate based 
5 on the received data according to equation (1). A weighted time window averager 
18 stores the output of the multiplier 16 and calculates a weighted average of the 
data-aided channel estimate according to equation (9). 

A channel estimate combiner 20 receives the output of the weighted time 
window averager 18 and the pilot symbol based channel estimate, calculates the 
10 variances of the DA channel estimate and the PA channel estimate, and 

generates the channel estimate according to equations (10), (13) and (14). The 
channel estimate is then used in the conventional manner to determine the 
transmitted symbols , and is also feedback to the LLR calculator 12 to be used 
in the LLR calculation. 

15 The invention being thus described, it will be obvious that the same may be 

varied in many ways. Such variations are not to be regarded as a departure from 
the spirit and scope of the invention, and all such modifications are intended to be 
included within the scope of the following claims. 
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